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p ■ We study the Higgs-boson mass spectrum of the minimal supersymmetric standard model, 

p ■ in which the tree-level CP invariance of the Higgs potential is broken explicitly by loop 

^ I effects of soft-CP-violating Yukawa interactions related to scalar quarks of the third genera- 

tion. The analysis is performed by considering the CP-non-invariant renormalization-group 
K> ! improved effective potential through next-to-leading order that includes leading logarithms 

due to two-loop Yukawa and QCD corrections. We find that the three neutral Higgs par- 
ticles predicted by the theory may strongly mix with one another, thereby significantly 
modifying their tree-level couplings to fermions and to the W^ and Z bosons. We analyze 
the phenomenological consequences of such a minimal supersymmetric scenario of explicit 
CP violation on the production rates of the lightest Higgs particle, and discuss strategies 
for its potential discovery at high-energy colliders. 
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1 Introduction 

Despite the great phenomenological success of the miniinal standard model (SM) at colhder 
and lower energies, its full experimental vindication is not yet complete. The Higgs boson 
iJ, the ultimate cornerstone of the SM responsible for endowing the observed fermions 
and the W^ and Z bosons with masses, still remains elusive thus far. Recent experiments 
at the CERN Large Electron Positron Collider operating at energies of 189 GeV (LEP2) 
place a severe lower bound on Mh, i.e. Mh > 95.5 GeV at 95% confidence level (CL) ||T|. 
On the other hand, global experimental data and theoretical analyses of radiative effects 
suggest that if nature indeed realizes the SM Higgs boson, it is then very unlikely that its 
mass be much larger than about 250 GeV 0. Nevertheless, there are many theoretical 
reasons to believe that the SM represents the low-energy limit of a more fundamental 
theory whose first clear signals are expected to be seen in experiments accessing energies in 
the range of 0.1 to 1 TeV. Especially, supersymmetry (SUSY) appears theoretically to be a 
compelling ingredient for a successful unification of gravity with all other fundamental forces 
in nature by means of supergravity and superstrings. In its minimal realization, the minimal 
supersymmetric standard model (MSSM), SUSY must be broken softly, in agreement with 
experimental observations 0. Unlike the SM, the MSSM offers an appealing solution to the 
gauge hierarchy problem, which is reflected by the perturbative stability of such a theory 
from the electroweak to the Planck energy scale. Because of the holomorphicity of the 
superpotential, the MSSM must contain at least two Higgs doublets, denoted as $i and 
$2, with opposite hypercharges, ^($2) = —^('^1) = 1, so as to give tree-level masses to 
both up and down families, and to cancel the triangle anomalies. 

Even though SUSY requires the presence of two Higgs doublets in the theory, the 
so-extended Higgs sector of the MSSM remains very predictive. This is because, at the tree 
level, all four-dimensional quartic couphngs in the Higgs potential are not independent, but 
related to the known electroweak coupling constants g^ and g' of the gauge groups SU(2)i 
and U(l)y, respectively. The CP-conserving MSSM predicts three Higgs states: two of the 
Higgs bosons, h and H, are even under CP and the Higgs boson A has CP-odd parity. 
Beyond the Born approximation, extensive theoretical studies based on renormalization- 
group (RG) methods and diagrammatic techniques have shown @,^J^ that the lightest 
CP-even Higgs boson, h, must possess a mass below 130 GeV. This upper bound is reached 
for large values of the ratio of Higgs vacuum expectation values, tan/? = ($2)/(*&i) > 15. 
For low values of tan/5 ~ 2, the upper bound on the mass of h decreases substantially, 
i.e. Mh ^110 GeV. For such low tan/5 scenarios, the current experimental limit on Mh is 
almost equal to the SM Higgs-mass bound for large values of Ma, i.e. Mh > 95 GeV at 95% 



CL, decreasing slightly for low values of Ma-, of the order of the weak scale. Therefore, the 
present experimental bounds put strong restrictions on models with low values of tan/3, 
close to the infrared fixed-point value 0J§]. Therefore, next-round experiments at LEP2 
turn out to be very crucial, as they can potentially exclude a significant portion of the 
parameter space of the CP-conserving MSSM |§] . 

Most analyses of the Higgs-boson mass spectrum of the MSSM have been performed, 
in the existing literature, only within the restricted framework of an effective CP-invariant 
Higgs potential. Recently, it has been shown |Q, however, that the tree-level CP invari- 
ance of the MSSM Higgs sector can be broken sizeably by one-loop effects that involve 
trilinear CP-violating couplings of the scalar top and bottom quarks to the Higgs states. 
As a consequence, the high degree of the tree-level mass degeneracy between H and A 
may be considerably lifted at one loop. Within the context of general two-Higgs-doublet 
models, the latter possibility has been extensively discussed by several authors, in con- 
nection with observable phenomena of resonant CP violation through HA mixing [|10| at 
future high-energy colliders, such as the Next Linear e^e~ Collider (NLC) |TT|,[T2|, the 
CERN Large Hadron Collider (LHC) ]T0|,|13[ and the proposed First Muon Collider (FMC) 



||1J,[T5[. Another important consequence of Higgs-sector CP violation in the MSSM is that 
the loop-induced CP-violating hA mixing can be of a size comparable to Mh, which may 
affect the predictions obtained for the mass of the lightest Higgs boson. CP violation and 
a light neutral Higgs boson are essential ingredients to account for the observed baryonic 
asymmetry in the Universe, through the mechanism of electroweak baryogenesis in the 
MSSM P . 



In this paper, we shall systematically study the mass spectrum of Higgs bosons in 
the MSSM with explicit CP violation. In such a scenario, both Higgs-boson masses and 
their couplings to fermions, W^ and Z bosons are significantly affected by the presence of 
CP-violating interactions. Therefore, we shall pay particular attention to the predictions 
obtained for the production rates of the lightest Higgs boson at LEP2 and other high-energy 
machines, such as the upgraded option of the Tevatron at Fermilab. Our analysis will be 
based on the computation of the CP-non-invariant RG-improved effective potential up to 
next-to-leading order. The dominant contributions to the effective potential come from the 
top (t) and bottom (6) quarks, as well as from their supersymmetric partners. For this 
purpose, we neglect chargino and neutralino quantum corrections. However, we include 
leading logarithms due to two-loop QCD and t and h Yukawa corrections ||T7| , [T8| . As has 
been explicitly shown in [0J§], these corrections improve the effective Higgs potential by 
minimizing its scale dependence significantly. 

CP-violating low-energy constraints, especially those coming from the electric dipole 



inoinent (EDM) of the neutron and the electron, play an important role in our analysis 
|T9| , [20| , |2ll , |22| , |23| , |24| , p5| . However, there have been several suggestions to evade these con- 
straints, without suppressing the CP- violating phases of the theory. One option is to make 
the first two generations of scalar fermions rather heavy, e.g. above 1 TeV |2^. Another 
possibility is to arrange for partial cancellations among the different EDM contributions 
||23| . Finally, it is possible to make the quantum corrections of the first two generations to 
EDMs negligible by requiring a kind of non-universality in the soft-trilinear Yukawa cou- 



plings [Q. Nevertheless, the aforementioned options do not prevent the supersymmetric 
analog of the two-loop Barr-Zee mechanism for generating EDMs from becoming large for 
high values of tan/? [^. Therefore, in our study, we shall consider these last 'direct' EDM 
constraints related to the CP- violating phases of scalar quarks of the third generation. 

The organization of the paper is as follows: in Section 2 we consider the CP- violating 
RG-improved effective Higgs potential of the MSSM, and derive the minimization conditions 
related to the Higgs ground state. In addition, we calculate the general mass matrices of 
the neutral and charged Higgs bosons. Technical details are relegated to the Appendix. In 
Section 3 we compute the three mass eigenvalues of the (3 x 3)-dimensional mass matrix of 
the neutral Higgs bosons, and the respective mixing angles. In Section 4 we consider the 
effect of EDM constraints on the CP-violating parameters related to the third-generation 
scalar fermions. Section 5 discusses the interactions of the Higgs particles with fermions, 
and with the W^ and Z bosons in the presence of CP violation. Furthermore, we analyze 
the phenomenological implications of these interactions for Higgs-boson searches at LEP2, 
and discuss the prospects of probing such a minimal SUSY scenario of explicit CP violation 
at an upgraded Tevatron machine. Section 6 summarizes our conclusions. 

2 The CP-violating Higgs potential of the MSSM 

The MSSM introduces several new parameters in the theory that are absent from the SM 
and could, in principle, possess CP-odd phases [|1^. Specifically, the new CP-odd phases 
may come from the following parameters: (i) the mass parameter fi, which involves the 
bilinear mixing of the two Higgs chiral superfields in the superpotential; (ii) the soft-SUSY- 
breaking gaugino masses mx, where A collectively denotes g, W and B, i.e. the gauginos 
of the gauge groups SU(3)c, SU(2)i and U(l)y, respectively; (iii) the soft bilinear Higgs- 
mixing mass wi^g; which is sometimes denoted as Bfi in the literature; and (iv) the soft 
trilinear Yukawa couplings Af of the Higgs particles to scalar fermions. If the universality 
condition is imposed on all gaugino masses at the unification scale Mx, then mx has a 
common phase. Likewise, the different trilinear couplings Af are all equal at Mx, i.e. 



Af = A. Here, one may slightly deviate from exact universality by assuming that A is 
a matrix in the flavour space |24]. In particular, it has been argued recently ^5\ that 



many dangerously large contributions to the electron and neutron EDMs may naturally be 
avoided by choosing the trilinear coupling of the Higgs fields to the scalar quarks of the 
first and second generation to be much smaller than the one of the third generation, i.e. 
A^^ (0,0,1) A 

It is known that the conformal-invariant part of the supersymmetric Lagrangian pos- 
sesses two global Peccei-Quinn-type symmetries: 

(i) The U(1)q symmetry, with Q assignments Q{Hi) = 1, Q{H2) = —2, Q{Q) = Q{L) = 
0, Q{U) = 2 and Q{D) = Q{E) = —1, where the caret on the fields symbolizes 
superfields. Hi and H2 are the Higgs superfields, which have opposite hypercharges 
Y{H2) = —Y{Hi) = 1, and Q (L), U and D (E) are the chiral multiplets of the 
quark (lepton) left-handed doublet, the right-handed up quark and the right-handed 
down quark (lepton), respectively. The chiral multiplets carry the hypercharges: 
Y{Q) = 1/3, Y{L) = -1, Y{U) = -4/3, Y{D) = 2/3, and Y{E) = 2. The U(1)q 
symmetry is broken by the fi and r/i^g parameters. 

(ii) The U{1)r symmetry acting on the Grassmann-valued coordinates 6 and 6, i.e. e^°'6 
and e~*"6'. So, the 6 coordinate of superspace carries charge 1. Moreover, all matter 
superfields carry charge 1 and all Higgs superfields carry charge 0. Under such a 
transformation, the gaugino fields carry charge 1. The superpotential carries charge 
2. Consequently, this symmetry is broken by the Majorana masses of the gauginos 
as well as by the scalar-fermion-Higgs trilinear couplings Af and the parameter /i. 

As a consequence, not all phases of the four complex parameters {/x, mfg, '"^a, A} turn 



out to be physical |20[, i.e. two phases may be removed by redefining the fields accordingly. 
Employing the global symmetries (i) and (ii), one of the Higgs doublets and the gaugino 
fields A can be rephased in a way such that mx and m^g become real numbers. As we will 
see, the fact that mfj is made real complies also with the CP-odd tadpole constraint on 
the Higgs potential at the tree level p. Thus, arg(/i) and aig^A) are the only physical 
CP-violating phases in the MSSM supplemented by universal boundary conditions at low 
energies. 

Denoting the scalar components of the Higgs superfields Hi and H2 by $1 = zr2$^ 
{t2 is the usual Pauli matrix) and $2, the most general CP-violating Higgs potential of the 

* Observe that, owing to the RG evolution, even starting with universal boundary conditions at high 
energies, the low-energy parameters tend to be non-universal. 



MSSM may conveniently be described by the effective Lagrangian 

+ Xs{<l>\<l>i){<l>l<l>2) + A4($I<l'2)($5$i) + Asl^I'l's)' + A*(<l>5$i)2 (2.1) 

+ A6($I$l)($l$2) +A*($I<l>i)(<l>5$i) + A7($^$2)($I$2) +A;($^$2)($^$i). 

At the tree level, the kinematic parameters are given by 

2 2||2 2 2||2 \ \ -'-/2|/2\ 

A3 = -\i9l-9''), A4 = ^gl, As = As = A7 = 0. (2.2) 

In Eqs. ( p.l|) and ( p.2|) , ?7i^, r7i2 ^i^d m^2 ^'^s soft-SUSY-breaking parameters related to 
the Higgs sector. Beyond the Born approximation, the quartic couplings A5, Ag, Ay receive 
significant radiative corrections from trilinear Yukawa couplings of the Higgs fields to scalar- 
top and scalar-bottom quarks. These parameters are in general complex. The analytic 
expressions of the quartic couplings are given in the appendix. 

Our next step is to determine the ground state of the MSSM Higgs potential. To this 
end, we consider the linear decompositions of the Higgs fields 

V 75 (^1 + '^1 + ^"1) / V 73 ^^2 + </>2 + ^02) J 

where f 1 and f 2 are the moduli of the vacuum expectation values (VEVs) of the Higgs dou- 
blets and ^ is their relative phase. Without any loss of generality, in the parameterization 
of the Higgs doublets in Eq. (|2.3|), we have adopted a weak basis in which the VEV f 1 {V2) 
and the quantum fluctuation 0i (^2) have the same phase. Furthermore, we assume the 



absence of any CP-odd component due to spontaneous CP violation [^]. Although radia 



five corrections can, in principle, lead to a spontaneous breakdown of CP invariance in the 



MSSM p^, such a particular scenario, however, predicts an unacceptably small mass for 



the CP-odd Higgs scalar A, i.e. Ma < 40 GeV, and it is therefore ruled out experimentally 



The VEVs vi and V2 and the phase ^ can now be determined by the minimization 
conditions on Cy This is achieved by requiring that the following tadpole parameters 
vanish: 

ar 1 

T^i = ("qT~) = ^c^{/"i + Re(mi2e*^) tan/? + t;^ Aic^ + - (A3 + A4)4 

+ Re(A5e2*«)4 + ^Re{\ee'^)spCp + ^ Re(A7e'«)4 tan/?] } , (2.4) 
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^sp{f4. + Re(mi2e*^) cot/3 + f^ A2S^ + - (A3 + A4)c; 



+ Re(A5e2*«)4 + ^Re{Ke'^)clcot P + ^ReiXre'^)spCf3] } 
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+ ^ Im(A7e*«)f 24 
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:2.6) 






- Im(A7e^«).;24 



(2.7) 



where s^ = sinx, c^^ = cos a;, tan/5 = ^2/^1 and v'^ = vf + v^. Furthermore, we assume the 
absence of charge-breaking minima, i.e. variations of £y with respect to (pf and 02^ vanish 
identically. It is now easy to see that the orthogonal rotation of the CP-odd fields. 




cos /3 — sin /3 
sin P cos P 




(2. 



gives rise to a flat direction of the Higgs potential with respect to the G° field, i.e. 
{dCv/dG^) = 0. In the newly defined weak basis, the G° field becomes the would-be 
Goldstone boson, which is absorbed by the longitudinal component of the Z boson. More- 
over, the orthogonal rotation (|2.8|) of the CP-odd fields yields a non-trivial CP-odd tadpole 
parameter: 

1 
2 

(2.9) 



T 



^ da ' 



— V 



Imiml^e'^) + ^(Ase^^^) v'^spCp + ^ lm{XQe'^y cl 



In the CP-invariant limit of the theory, both Ta and the phase ^ vanish. Since m\2 is taken 
to be real at the tree level, a non-zero value of the phase ^ is first generated at the one- 
loop level 1^. Nevertheless, in a general two-Higgs-doublet model with Higgs-sector CP 
violation, the phase ^ already occurs in the Born approximation. For the sake of generality, 
we shall keep the full ^ dependence in the analytic results. 

It is now interesting to discuss the conditions under which the Higgs sector of the 
MSSM respects the CP symmetry. We find that the CP invariance of the Higgs potential 
is assured only if 



Im(mi2A5) = Im(mi2A6) = Im(mi2A;) = 



(2.10) 



If we assume a kind of universality between the trilinear Yukawa couplings at low energies, 
A = At = Ab, and neglect the small chargino and neutralino contributions, the phases of 



the quartic couplings A5, Ae and A7 are then related to one another. Employing the analytic 
results of the quartic couplings given in the Appendix, it is easy to show that 

Im(A*A^) = Im(A*A?) = 0. (2.11) 

However, even in this case, the phase of the complex soft parameter m^2 is not restricted 
by any universal boundary condition imposed by minimal supergravity models at the uni- 
fication point. In other words, the CP invariance of the Higgs potential of the MSSM holds 
true only if the condition 

lm{mllfiA) = (2.12) 

is satisfied. Within the most general framework of the MSSM, the equality ( |2.12| ) [or 
equivalently ( p.lO| )] can be violated, thus giving rise to observable CP violation. 

CP violation in the Higgs potential of the MSSM leads to mixing mass terms between 
the CP-even and CP-odd Higgs fields [§]. Thus, one has to consider a (4 x 4)-dimensional 
mass matrix for the neutral Higgs bosons. In the weak basis (G°, a, 0i, ^2), the neutral 
Higgs-boson mass matrix A^q may be cast into the form 



Ml M] 
Mlp Mi 



Ml = I 7J Z7 I . (2-13) 



where Ml and M^ describe the CP-conserving transitions (G"^, a) -^ (G^, a) and 
(01)02) —* (01,02), respectively, and Mlg = (M^gp)^ contains the CP-violating mixings 
{G^, a) ^^ (01, 02)- The analytic form of the sub matrices is given by 

V V 

spT^^ - CfjT^^ , ,2 S/3 tan (5 T^, + cp cot (3 T^ 



■^P - I .^^r. -r.T. _ s p tsin p\, + Cf, cot (3 T^, I ' ^^.14) 

V a V 




A^2 _ ,2 1 Im(A5e2^«)s/3 + Im(A6e^«)c/3 , ^a . ^p ^p . .. , .^ 



(2.16) 

_ 2 ; 2Ai4 + 2Re(A5e2*«)4 + 2Re(A6e^«)s^c^ X-uSpCfs + Re{X6e'^)c} + Re(A7e^«)s^ 
A34S/3C^ + Re(A6e^«)4 + Re(A7e^«)4 2A24 + ^M^^e^'^H + 2Re{X7e'^)spCp 

In Eqs. (|2.14|) and ( p.l6D , we have used the abbreviations A34 = A3 + A4 and 




Mf 



— {Re(m?2e'«) + v^ [2Re(A5e''«)s;3C^ + ^ Re(A6e^^)4 + ^ Re(A7e'«)4] } . (2.17) 



S/3C/3 



If the MSSM is invariant under the CP symmetry, Ma is then the physical mass of the 
CP-odd Higgs scalar [§. 

Correspondingly, the charged Higgs-boson mass matrix }A\, spanned in the mass 
basis (G^, H^), may be obtained by the Lagrangian 






{G+, H+) M 



G- 
H 



(2.18) 



with 



Ml 



V 
V ~^ V 



SjjT^i — cpT^. 



•2 -T, 

- — I 



V 



<± 



S/3 tan /5 T(^^ + C/3 cot (3 T^^ 

V 



(2.19) 



and 



M^ 



1 



H± 



{Re(7 



SpCp 

+ -Re(A7e^« 



mfac'^) + v" 



- \ASpCp + Re{\^e^'^)s(iCfi + - Re(A6e*^)c; 



)sl]} 



(2.20) 



From Eqs. ( p. 17] ) and (|2.2(]| ), we observe that M^± is related to the mass of the would-be 
CP-odd Higgs scalar Ml through 



Mt 



Ml± - -X^v^ + Re{\r,e^'^y 



(2.21) 



Taking this very last relation into account, we may express the neutral Higgs-boson masses 
as functions of Mh+, /i, A, ^fe, the common SUSY scale Msusy, tan/3 and the physi- 
cal phase ^. Since we neglect chargino and neutralino contributions, we can absorb the 
radiatively induced phase ^ into the definition of the /i parameter. 

It is worth stressing that, even though CP violation decouples from the sector of the 
lightest Higgs boson for large values of M11+ ~ M^, this decoupling property of CP non- 
conservation does not generally persist for the system of the two heaviest Higgs bosons. 
This point may formally be seen as follows: in the large M^ limit, assuming that the quartic 
couplings are kept fixed, the submatrix A^| given in Eq. ( p. 161) has one mass eigenvalue 
which approaches M^, and corresponds to the heaviest CP-even Higgs boson, while the 
other one is small at the electroweak scale, related to the lightest CP-even Higgs boson. 
Furthermore, the submatrix M.'\, in Eq. (|2.14| ) has only one non-zero mass eigenvalue equal 
to M^, corresponding to the mass of the would-be CP-odd Higgs scalar. Thus, for large 
Mf values, it is easy to see that the effective (2 x 2)-dimensional submatrix, which may 
be formed by the would-be CP-odd and the heaviest CP-even Higgs bosons, also contains 



off-diagonal CP- violating terms coming from Mgp in Eq. ( p.l5| ). These off-diagonal CP- 
violating matrix elements are generically of the order of the difference of the diagonal entries 
of the effective submatrix, thereby giving rise to a strong mixed system of CP violation. 
As we will see in Section 5, numerical estimates of CP violation in the heavy Higgs sector 
offer firm support of this observation. 



3 Higgs-boson masses and mixing angles 

In this section, we shall evaluate the physical masses of the neutral Higgs bosons and the 
mixing angles related to the diagonalization of the general 4x4 matrix TMq in Eq. (|2.13|) . 
Even though our primary interest is in the CP-violating MSSM, the validity of the analytic 
expressions that we shall derive here extends to the most general class of CP-violating two- 
Higgs-doublet models. After setting all tadpole parameters to zero, we easily see that G° 
does not mix with the other neutral fields and so becomes an independent massless field, 
as it should be on account of the Goldstone theorem 0. Then, A^q effectively reduces to 
a (3 X 3)-dimensional matrix, Ai%, which is spanned in the weak basis (a,0i,02)- 

The mass eigenvalues of A^^ are obtained by solving the characteristic equation of 
cubic order 

x^ + rx^ + sx + t = 0, (3.1) 

with 



MM%), 



Tt\M 



N) 



Tt{M 



N^ 



-det{M%). 

To this end, it proves useful to define the following auxiliary parameters: 

3s - r^ 



+ t. 



(3.2) 



p = 


3 




2^3 ^ 


q = 


27 


D = 


p3 2 

27 4 



(3.3) 



To ensure that the three eigenvalues are positive, it is necessary and sufficient to require 
that 

D < 0, r < 0, s > 0, t < 0. (3.4) 
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Imposing these inequalities on the kinematic parameters of the theory, we may express the 
three mass eigenvalues oi J^j^ as 

pI = -7;r + 2 \/-p/3 cos ( ^ 

pI 
pI ^ 

with 



3 ■ V -' V 3 y ' 

-^r + 2y^cos(| + y), 

-^r + 2y^cos(| - y), (3.5) 



(f = arccos ( ) . (3.6) 

^ 2J-py27^ 



Since the Higgs-boson mass matrix Ai% is symmetric, we can diagonalize it by means 
of an orthogonal rotation O as follows: 

0^M%0 = diag(Mi3, Ml^, M^J . (3.7) 

Note that some arbitrariness exists in assigning the mass eigenvalues pi, p2 and ps in Eq. 
( p.5|) to those related to the mass eigenfields Hi, H2 and H3 in Eq. ( p.7|) . For clarity of the 
presentation, we define these fields such that 

Mh, <Mh, <Mh,. (3.8) 

Alternatively, these fields could be defined in such a way that we have, in the CP-invariant 
limit of the theory. Hi = h, H2 = H, H^ = A, where h and H denote the lightest and 
heaviest CP-even Higgs bosons, respectively, and A is the CP-odd Higgs scalar. However, 
we should use the former definition, as the latter often leads to discontinuities in the Hi 
mass values, when plotted as a function of the MSSM parameters. 



In general, the orthogonal matrix O in Eq. ( |3.7| ) can be described in terms of the 
three physical Euler-type angles x? V' and 6. We parameterize O, assuming x,?/) — * or 
7r/2 in the CP-conserving limit of the theory, as follows: 



O 



/ c^c^ 


-S^Cg - S^C^Sg 


Sx-^e - S^C^Cg \ 


s^c^ 


c^ce - s^s^se 


-C^S0 - S^S^Cg 


\ H 


SgC^ 


cec^ J 



(3.9) 



It is convenient to find first the entries Oij, and then determine the three rotational angles 
by the obvious relations: 

ijj = arcsin (O31) , x = arcsin ( -), 6 = arcsin ( -), (3.10) 

^ cos 1p ^ ^ COS 'Ip ' 
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where the mixing angles take on values in the interval (— 7r/2, 7r/2]. 

If M^p with z, j = 1, 2, 3, denote the matrix elements of Aij^, the elements Oij can then 
be obtained by appropriately solving the underdetermined coupled system of equations, 

EkMlO,, = Ml^^^^p,f. 



(Mf, - Mfj^^JOu + Mt^02^ + Mt,Os^ 

33 



MlOu + MI2O2. + (M. 
More explicitly, we have 



M' 



H(4_,))C'3i 



0, 
0, 
0. 



(3.11) 
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( la^il/Ai 
2/i/Ai 
V ^i/Ai 



a;2/A2 X3/A3 \ 
1^/2 1 /A2 ys/As 

2:2/ A2 izsi/As ; 



(3.12) 



where 
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Mh 



12 
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Mh 
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Mh 
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\V2\ 
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Ai = ^Xl + 1/2 + z2 
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' Mh 
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'Xl 



(3.13) 
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Mh - Ml, 


' Mh 
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-Ml, Mh 


Mh 


Mh-Mh 



(3.14) 



In Eq. ( p.l4| ), we have used the abbreviation s^^ = sign(x), which is an operation that 
simply gives the sign of a real expression x. Notice that the parameterization of O in terms 
of Xi, yi, Zi may be chosen, in a way such that indefinite expressions do not occur in the 
CP-conserving limit of the theory. 



4 EDM constraints 

As we mentioned in the introduction, the EDM of the electron and the neutron pro- 
vide the most stringent constraints on the CP-violating parameters of the MSSM 
19| , |20| , pll , ^ , ^^ , p5| . There have been several suggestions to suppress the EDM contribu- 



tions coming from the first two families of scalar quarks without making the CP-violating 
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Figure 1: Two-loop contribution to EDM and CEDM of a light fermion / in the CP- 
violating MSSM (mirror-symmetric graphs are not shown). Note that f does not contribute 
to the CEDM of a coloured fermion /. 



phases of the theory very small p2| , p3tp4| , p5| . To be specific, the following three possibilities 
may be considered: one can make the first two families of scalar fermions rather heavy, 



having a mass of order few TeV [^ . Another option is to arrange for partial cancellations 
among the different EDM contributions p3| . In this case, the CP- violating phases of the 
theory turn out to be rather correlated. Finally, an interesting alternative is to adopt a 



slightly non-universal scenario for the trilinear couplings Aj p^ , |25[] . In particular, one may 
require that arg(/i) < 10~^ and Af = (0,0, 1)^4 [Q. In the latter scheme, A^- = At = Ah 
are the only large trilinear couplings in the theory with CP-violating phases of order unity. 



Furthermore, assuming that gluinos are heavier than about 400 GeV p3[, one may signifi- 
cantly reduce the size of the EDM effect due to Weinberg's three- gluon operator 0] well 
below the present experimental bound. 

In the aforementioned SUSY scenarios, however, the two-loop Barr-Zee-type [^ 



contribution to the electron and neutron EDMs, shown in Fig. |1], can still be potentially 



large ||2^ . It is important to notice that the very same kind of loop graphs generated by 
scalar top and bottom quarks are also responsible for CP violation in the Higgs sector 
of the MSSM. Therefore, we are compelled to take these two-loop EDM constraints into 
consideration. The graphs displayed in Fig. |I| give rise to both an EDM {dj/e) and a 
chromo-EDM (CEDM) [dfi / 9s) of a light (coloured) fermion /, i.e. 
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(4.2) 



where Qj (Qf) stands for the electric charge of a (scalar) fermion given in |e| units, Nj is 



the colour factor of the scalar fermion / [N? = 1 and A^^ 



m 



3), Rf = cot/3 for the 



up- family fermions, Rf = tan/5 for the down-family ones. In addition, we have defined 

sin 29fmflm{fie^^f ) 



Fiz) 



Rf 



dx — 
z 



sin P cos P v^ 
x{\ — x\ 



x{\ 



In 



X] 



x(\ — x) 



(4.3) 
(4.4) 



In Eq. ( ^.3| ), 5/ = arg(y4j — RffJ,*) and 6f indicates the mixing angle between weak and 
mass eigenstates of /. Since the off-diagonal elements of the scalar-quark and lepton mass 
matrices of the third generation may be larger than the difference of the diagonal entries, 
angles 6f close to 45° are obtained in a natural way. Further discussion and more details 
of the calculation may be found in [^ . 



The present status of measurements of the electron and neutron EDMs is as follows 



m 



/4 

^ e 

(7) 

/ dn 

^ e 



exp 



exp 



exp 



(-0.27 ±0.83) X 10"^^ cm [33], 
(0.18 ± 0.12 ± 0.10) X 10"^^ cm [§3 
(0.26 ± 0.40 ± 0.16) X 10"^^ cm [|5|. 



(4.5) 
(4.6) 
(4.7) 



The experimental numbers listed here contain an amount of theoretical uncertainty origi- 
nating from the model used for the heavy atoms, such as ^°^T1, or from the description of 



neutron's wave function in the heavy nucleus [^. Notwithstanding the possible uncertain- 
ties in the determination of d^, we should regard the la upper bound, \de/e\ < 1.1 x 10~^^ 
cm, stated in Eq. ( [4.5|) as a conservative one, when compared to the improved bound in 
Eq. ([4.6|) . In particular, the 2a upper bound on the electron EDM coming from the latter 
experimental analysis is \de/e\ < 0.5 x 10~^^ cm. Finally, the Icr and 2a upper bounds on 
the neutron EDM are |(i„/e| < 0.69 x 10"^^ cm and \dn/e\ < 1.12 x 10~^^ cm, respectively. 

It is interesting to confront the theoretical predictions obtained for d^/e and dn/e in 
the MSSM with the corresponding experimental bounds mentioned above. In Table |l|, we 
present numerical estimates for the case \A\ = \At\ = \Ab\ = \Ar\ = 1 TeV, arg(A) = 90° 
and Ma = 150 GeV, while the kinematic parameters tan/5, /i, Msusy have been varied 
discretely. We should note that the largest contribution to the neutron EDM comes from 
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Table 1: Numerical predictions for the electron and neutron EDMs in the MSSM, using 
\A\ = \At\ = \Ab\ = \A^\ = 1 TeV, arg(A) = 90° and M^ = 150 GeV {mb{Mz) = 3 GeV). 

the CEDM of the d quark. Furthermore, scalar-top quarks have the biggest quantum effect 
on the e and n EDMs for 2 < tan/5 < 30. Taking the numerical estimates in Table ^ into 
account, we shall consider the following three representative scenarios: 



I. 


^SUSY 


= 0.5 TeV , 


A 


II. 


^SUSY 


= 0.5 TeV , 


A 


III. 


^SUSY 


= 0.5 TeV, 


A 



ITeV, /i = 2TeV, tan/5 = 2, 
ITeV, /i = 2TeV, tan/5 = 4, 
ITeV, /i = lTeV, tan/3 = 20, 



(4.^ 



In the next section, we shall analyze the phenomenological consequences of the scenarios 
given by Eq. ( [4.8|) on the currently operating and future high-energy colliders. 

5 Phenomenological implications for Higgs searches 

We shall discuss the main phenomenological consequences of CP violation in the Higgs sec- 
tor of the MSSM on the Higgs-boson mass spectrum and on the production cross sections 
of the lightest Higgs boson Hi. In general, CP violation modifies the couplings of the Higgs 
particles to fermions and to the W and Z bosons, as well as their self-interactions. Further- 
more, quadrilinear interactions involving Higgs bosons change as well, when CP-violating 
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Figure 2: Numerical estimates of (a) Mh-^ and (b) 2\Mh2— Mh-j,\/ {Mh2 + ^Hz) as a function 
of the CP- violating phase arg(At). 
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effects due to a mixing of Higgs states are included. However, the latter interactions as well 
as the trilinear Higgs self-couplings are generally sub-dominant in production processes, and 
hence we should not consider these here. 

At the high-energy machines LEP2 and Tevatron, the dominant production mecha- 
nism of the lightest Higgs boson is the Bjorken process in association with W and Z bosons 
|37| , |38| , e.g. e"'"e~ -^ HiZ or the partonic process mJ — > W^Hi. Such reactions involve the 
couplings of Higgs bosons to W and Z bosons. In the presence of CP violation, these 
couplings may be read off by the Lagrangians 

Chvv = 9n.Mw{cfs02^ + SfsOs,)(H^^_,)W;^W-'^ + ^H^^_,)Z^Z^), (5.1) 
Chh±w^ = y(c/303. - Sf,02^ + tOu)W+'^{H^^.i^td^H-) + H.c, (5.2) 



J-HHZ 



On (cfsOsj - sp02j) - Oij (cfsOsi - spO 



2i) 



4c„ 

xZ'^(if(4_.)9^i^(4-,)), (5.3) 

■^-> — > < — 

where Cy^ = Mw/Mz and d^ = d^ — d^. Note that the Z boson can only couple to 
two different Higgs particles as stated in the Lagrangian ( p.3|) . The reason is that Bose 
symmetry forbids any antisymmetric derivative coupling of a vector particle to two identical 
real scalar fields. 

Making now use of the following identity, which governs the matrix elements of O 

(assuming detO = 1): 

1 3 

n,m,i,j=l 

it is not difficult to derive an important relation between the couplings of the neutral Higgs 
bosons to the gauge bosons, namely 

gHkVV = ^ijkQHiHjZ, (5.5) 

where gUiVV iY = W"^, Z) are the Higgs-gauge-boson couplings normalized to the SM 
value, i.e. 

gH.vv = cp02i + spO^i, (5.6) 

while gUiHiZ is defined by the expression between the brackets in Eq. (|5.3|), 



gH,H,Z = Oii{cpO'ij - Sp02j) - Oij{cpO-ii - Sp02i). (5.7) 

From the relation (|5.5| ) and the unitarity constraint 

Hal^zz = 1, (5.8) 



j=i 
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the immediate result is that the knowledge of two gu^zz is sufficient to determine the whole 



set of couplings of the neutral Higgs to the gauge bosons |^ . 

The essential difference between a CP-conserving MSSM Higgs sector and a CP- 
violating one is that mixing effects between the would-be CP-odd and CP-even Higgs 
bosons are present in the latter case. Such scalar-pseudoscalar mixing effects are induced 
by radiative corrections to the Higgs potential |^. The characteristic size of the CP- 
violating off-diagonal terms in the Higgs-boson mass matrix may be estimated by 

|2 I, ,12 



x|6..^^. >...r.'. . :t:»." i. (5^9) 



\At\ \[i\ sin20cp l/^M 

M|usY ' tan/?M|usY' sin0cp ^susy 
where the last bracket summarizes the relative size of the different contributions, and 



<\>cv = arg(A/i) + e. (5.10) 

For |/i| and \At\ values larger than the arithmetic average of the scalar-top-quark masses 
squared, denoted as Msusy (cf. Eq. ( |A.l(j| )), the CP-violating effects can be significant. 
For instance, if \fi\ ^ \At\ ^ 2Msusy, and 0cp — 90°, the off-diagonal terms of the neutral 
Higgs-boson mass matrix may be of the order of (100 GeV)^. These potentially large mixing 
effects have important consequences, since they lead to drastic variations in the definition 
of the neutral Higgs-boson masses and in the couplings of the Higgs states to the gauge 
bosons. Because of the quantum nature of the CP-violating mixing effects and the known 
decoupling property of heavy states in the loop in SUSY theories, the phenomenology of 
the lightest Higgs boson will only be important for low values of Mh+ ■ For the same reason, 
Mh+ values much larger than the electroweak scale lead to predictions for the mass of the 
lightest Higgs boson Hi and for the couplings of the Hi scalar to the gauge bosons which 
are equivalent to those obtained in the CP-invariant theory. The only difference in the 
CP-violating case is that the relevant scalar-top mixing parameter entering the definition 
of Mhi is now given by 

\At\ = \At-fj.*/tanP\. (5.11) 

In addition to the effects induced by scalar-top quarks, Yukawa interactions due to 
scalar-bottom quarks can also be significant for large values of tan/5, so as to lead to 
sizeable contributions to the elements of the Higgs-boson mass-matrix. As was discussed in 
Section 4, however, the contributions of the scalar-bottom sector are limited by constraints 
that originate from the electron and neutron EDMs. In general, unless cancellations occur 
between different contributions to the EDMs, the CP-violating quantum effects coming 
from the scalar-bottom sector are expected to be small. 
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Figure 6: The same as in Fig. 0, but with tan/3 = 4. 
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In Fig. 1^, we present numerical estimates of the lightest Higgs-boson mass Mh^ and 
of the relative mass splitting of the two heaviest Higgs bosons, H2 and H^, as a. function 
of (pcv = arg(^t)5 fo^ different values of the charged Higgs-boson mass and for tan/9 = 2, 
^susY = 0.5 TeV, \At\ = 1 TeV and /i = 2 TeV. We must remark that relatively large values 
of the \fi\ and \At\ parameters have been chosen, in order to make the CP- violating effects 
relevant. In fact, for these values of the trilinear couplings, we are at the edge of the limit of 
validity of the expansion we used to compute the expression of the quartic couplings in the 
RG-improved Higgs potential (see appendix). However, we may verify that the neglected 
terms would lead to a slight increase of the CP-violation effects we have computed. For the 
region of parameters we are considering, the omitted terms are of the order of finite two- 
loop corrections to the effective potential, which are not included in the analytic expression 
of the one-loop RG-improved effective potential. Moreover, if all SUSY-mass parameters 
are rescaled by a factor 2, the validity of the expansion improves dramatically, while the 
physical results are only slightly modified. Later on, we shall explicitly demonstrate this 
point by giving a specific example (cf. Fig. P). In this context, we believe that our results 
should be regarded as conservative estimates of the possible size of the CP- violating effects 
in the Higgs sector. 

In Fig. I, the dependence of Mh^ on aig^At) may be understood as follows: if the 
moduli of the parameters At and fi are kept fixed, the mass of the lightest neutral Higgs 
boson Hi starts increasing, as a function of a.Tg{At), up to a maximum; it then decreases 
rapidly. This kinematic dependence may be attributed to the fact that the modulus of the 
scalar-top mixing parameter At increases monotonically for the range < arg(Ai) < 180°. 
We find that the behaviour of Mh^ as a function of aTg{At) is very analogous to the one 
that would have been obtained for the lightest CP-even Higgs boson in the CP-invariant 
theory if we had varied \At\ given in Eq. ( p.ll| ). However, unlike the CP-conserving case. 



the effects of the Higgs-sector CP violation can give rise to much larger mass splittings 
between the two heaviest Higgs bosons H2 and H3 as a function of arg(y4i) [^. This fact 
is difficult to infer from the kinematic dependence of Mh^ alone. In Fig. |^(b), we see that 
even in the case of Mjj+ = 200 GeV, the CP-violating effects may lead to a relative mass 
splitting of the two heaviest neutral Higgs bosons of the order of 30%. 

In order to get an idea of how much of the above effects are due to the presence 
of CP violation, we plot in Fig. |^ the behaviour of the same neutral Higgs-boson masses 
as a function of the parameter /i, where 0cp = is considered and the scalar-top mixing 
parameter \At\ is fixed to the value that yields the maximum for the lightest CP-even Higgs- 
boson mass Mfi for large values of Mh+ . Figure ^ exhibits the dependence of M^ = Mh^ as 
a function of /x. We find that large values of /x tend to reduce Mh and increase the degree 
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of the mass splitting between the heaviest CP-even Higgs boson and the CP-odd Higgs 
scalar, 2\Mh2 ~ ^H3\/{Mh2 + ^Hs)- In particular, we observe that for the same value of 
Mhii 2\Mh2 — Mh-J/{Mh2 + Mhs) is much smaller in the CP-invariant theory than in the 
CP-violating one. 

It is worth investigating whether the CP-violating Higgs effects could lead to not-yet- 
explored open windows in the parameter space of the theory that cannot be accessed by 
the running experiments and are not present in the CP-conserving case. To this end, we 
plot in Figs. § and ^ the Higgs-gauge-boson couplings as a function of arg(yli), for the same 
choice of SUSY parameters as those in Fig. §. For low values of Mh+, the Higgs boson 
that couples predominantly to the Z boson is Hi {H2). For Mh+ = 140 GeV, there is an 
interesting region of parameters for which the HiZZ coupling becomes small, rendering the 
detection of the Hi boson via the HiZZ channel at LEP2 difficult. However, in the same 
region of parameters, the mass of the next-to-lightest neutral Higgs state, H2, is smaller 
than 95 GeV and the H2ZZ coupling becomes large. Therefore, in this particular kinematic 
range of parameters, the H2ZZ channel is the only relevant one that helps to exclude these 
small values of tan/5 and charged Higgs-boson masses at LEP2. For larger values of Mh+, 
the situation resembles more the CP-invariant theory. Even for Mh+ = 170 GeV, the 
suppression of the HiZZ coupling is not sufficiently large to produce a significant change 
in the present LEP2 bound. Therefore, very analogously to the CP-conserving case P], only 
a small region of parameters, for which the lightest CP-even Higgs boson is heavy enough, 
is still experimentally allowed, i.e. for large values of the scalar-top mixing mass parameter 
I Ail (arg(A() ~ 60°-90°) and for relatively large values of the charged Higgs-boson mass. 

More interesting is the situation at intermediate values of tan/3. Even though, as 
displayed in Figs. | and |^, the dependence of the Higgs masses on the CP-violating phases 
and the scalar-top mixing parameters is similar to the one obtained for tan/3 = 2, the region 
for which the coupling of the lightest neutral Higgs boson to the gauge bosons becomes small 
now displays larger values of the Hi and H2 masses. In this range of kinematic parameters, 
the experimental detection of the Hi and H2 particles becomes more difficult at LEP2. 
In Table |], we present the current reach of LEP2 for the i^i-boson mass as a function of 
QhiZZ, and independently, for {Mh^,Mh^) versus gHiH2Z, assuming Mh-^ ^ Mh2- As can 
then be seen from Fig. § for the case of a relatively light charged Higgs scalar Mh+ = 140 
GeV, we get regions for which the lightest Higgs-boson mass Mh^ is as small as 60-70 
GeV and the HiZZ coupling is small enough for the Hi boson to escape detection at the 
latest LEP2 run, with ^/s = 189 GeV. Moreover, the H2 boson is too heavy to be detected 
through the H2ZZ channel. In addition, either the coupling H1H2Z is too small or H2 
is too heavy to allow the Higgs detection in the H1H2Z channel (see also Fig. 110). To 
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9hiZZ 


Mh, 


9hiH2Z 


Mh, ^ Mh, 




[GeV] 




[GeV] 


1 


97 


1 


86 


0.6 


96 


0.9 


85 


0.5 


95 


0.8 


84 


0.4 


93 


0.7 


83 


0.3 


92 


0.6 


82 


0.2 


88 


0.5 


80 


0.1 


77 


0.4 


77 


0.08 


70 


0.3 


72 



Table 2: Present experimental sensitivity of the couplings gfj^zz ^"^^ 9hiH2Z ^s a function 
of Mhi , assuming Mh^ ~ M^z ■ The results are obtained from four experiments combined 
at LEP2 with ^/i = 189 GeV. 



better gauge the validity of the expansion that has been used for the quartic couplings in 
the effective Higgs potential (see also the discussion in the appendix), we have plotted in 
Fig. ^the same tan/5 scenario but rescaling Msusy, |A| = |^fe| and /i by a factor of 2, i.e. 
MsusY = 1 TeV, \At\ = |A| = 2 TeV and /i = 4 TeV. We then find that Mh^ has only 
slightly increased and the kinematic dependence of the HiZZ coupling on the CP- violating 
phase aig{At) remained almost unchanged when compared to Fig. p. 

For very large values of tan/3, the CP- violating effects in the Higgs sector are con- 
strained by the bounds on the EDM of the electron and neutron discussed in Section 4. If no 
cancellation mechanism is assumed to occur between the different EDM contributions, the 
CP- violating effects in the Higgs sector are then relatively small, and the phenomenological 
properties of the Higgs bosons become very similar to those obtained in the CP-invariant 



case. This feature is shown in Figs. |TT]-0, for the same quantities as those considered in 
Figs. ^|-^ Notice that the actual bound on tan/5 coming from EDM constraints strongly 
depends on the exact value of the CP-odd phase, which we choose to be equal to 90°. In the 
high-tan /5 regime, the running 6-quark mass rrih plays a central role. As we will see later 
on, nib turns out to be a model-dependent quantity in the MSSM. So as to make definite 
predictions, however, we have set mbijfit) = 3 GeV, where Wt is the top-quark pole mass 
(cf. Eq. (^;T|)). 

It is now interesting to discuss whether the Tevatron collider together with its future 
upgraded facilities have the potential capabilities of exploring the open windows that are 
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Figure 10: The same as in Fig. |^, but with tan/? = 4. 
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not accessible at LEP2 pO| , ^ . The Tevatron reach depends very strongly on the final 
run II luminosity. Furthermore, the reach is very much affected by a suppression of the 
Higgs-boson production rate in the HiVV channel. Hence, it is unlikely that the Tevatron 
will be able to observe the lightest Higgs boson in the regions where the effective HiVV 
coupling is suppressed. The search for the H2 boson looks more promising. Indeed, for 
high luminosities available at a later stage of the collider, e.g. 30 fb~^ per experiment, the 
Tevatron-discovery reach of the Higgs-boson mass can be as high as 125 GeV for a SM 
production rate and 115 GeV if there is a 0.7 suppression factor.[| From the results of 
Figs. |6|-|T0|, we can see that, as happens in the CP-conserving case [^], a high-luminosity 



Tevatron collider may be capable of covering most of the windows left open at LEP2. 

In Fig. |T5|, we analyze the variation of the coupling of the charged Higgs bosons H"^ to 



the lightest neutral Higgs boson Hi and the W"^ gauge bosons, for the same parameters as 
those considered in Fig. |^. The HiH^W~ coupling may be defined from Eq. ( [5.2| ) without 
including the weak gauge-coupling factor (7^/2, i.e. 

9h,h+w- = C13O33 - Sfs023 + iOis . (5.12) 

The HiH~^W~ coupling is relatively small when the coupling HiZZ is large, while it is 
enhanced when the coupling HiZZ is suppressed. On the other hand, a measure of CP 
violation in the HiH^W~ vertex may be obtained by analyzing the CP-odd quantity 
\hii{gjj^^+y^^)\/\gu^H+w~\'^ (see also Fig. |15|(b)). This CP-odd quantity shows a very in- 



teresting kinematic behaviour, as it can become of order 1 in large regions of the parameter 
space with Mh+ ~ 300 GeV. A consequence of large CP violation in the HiH^W~ coupling 
is that the decay rates for H^ —>■ HiW~^ and H^ —>■ HiW~ may become very different ^ . 



In the following, we consider the interactions of the neutral Higgs fields with the 
fermions. These interactions may be obtained by the Lagrangian 

3 

^Hf-f = - Y.H{A-i) \7ZTr^d{02i - ispOia^)d 

+ 7;7T^^'^{0u - icpOia^)u . (5.13) 

Obviously, the Higgs-fermion-fermion couplings are significant for the third-generation 
quarks, t and h. From Eq. (|5.13| ), we readily see that the effect of CP- violating Higgs 
mixing is to induce a simultaneous coupling of ifj, with z = 1, 2, 3, to CP-even and CP-odd 



fermionic bilinears [H^, e.g. to uu and ui'j^u. This can lead to sizeable phenomena of CP 



T These values rely on a combination of the HiWW and HiZZ channel for both experiments ]38| , 
assuming that the total branching ratio into b quarks remains almost unchanged with respect to the SM 
value. 
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and (b) 
'^\{9h2uu) {9H2uu)\/[{9H2uuf + {9H2uuf] as a function of arg(At). The definition of the cou- 
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plings Qii^dd and g^j^^^, with i = 1, 2, is given in Eq. (^T§. 
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violation in high-energy processes that involve decays of Higgs bosons into longitudinally 



polarized top-quark pairs |44 



For the discussion that follows, it proves convenient to define the following parameters: 

ai.uu = 033/Sf3 , gHruu = ^13 COt /? , g^H2Uu = 0^2/ Sf3 , Qh^uu = ^12 COt p , 

9H^dd = O23/C/3 , gH.dd = Oi3 tan/3 , gfj^^^ = O22/C/3 , Qn^dd = Ou tan/5 . (5.14) 

These parameters represent the scalar and pseudoscalar couplings of the Higgs bosons Hi 
and H2 to the up- and down-type fermions, normalized to the SM values. Then, the partial 
decay widths of Hi and H2 in the MSSM may be obtained by the SM ones, after multiplying 
the latter by the effective coupling factors [(fi-^i//)^ + (S'l/i//)^] and [igfj^ff)'^ + (fl-Ha//)^] 



(with f = u,d), respectively. In Figs. |T6| and [T^, we plot the effective coupling factors as 
a function of a.Tg{At). The behaviour observed is similar to the CP-conserving case. To 
be specific, the effective coupling factor related to the Hiff coupling (with i = 1,2) is 
close to the SM value, whenever the respective HiVV coupling approaches 1, while in the 
regions where the HiVV coupling is suppressed, the Hi coupling to down (up) fermions is 
enhanced (suppressed) by a factor tan/5 (1/tan/?). 

One may now construct quantities that can provide a realistic measure of CP violation 
in the HJf and ifa// couplings, i.e. 2\{gfj^jj) (^H,//)|/[(^i,//)^ + (^|^,//)^]. Figures |l| and 



|T9| exhibit the dependence of these CP- violating quantities related to the Hi and H2 bosons, 
respectively, as a function of aig{At). Figure ^ reveals that the CP- violating component of 
the Hiff coupling is large only for relatively light charged Higgs-boson masses, Mh+ < 180 



GeV. On the other hand, as can be seen from Fig. |T9|, the CP-violating component of the 
H2f f coupling may become of order 1 for a wide range of the parameter space, even for 
heavier charged Higgs-boson masses, e.g. Mu+ ^ 500 GeV. 

Apart from the CP-violating effects generated by the radiative mixing of the Higgs 
fields which we consider here in detail, there may be important CP-violating effects induced 
by one-loop vertex corrections. In the leptonic sector, these corrections are generally small 

Because of the large Yukawa and colour-enhanced QCD interactions, however, the 



radiative corrections to the tree-level 6-quark couplings to the Higgs bosons [^ may be 
important when the relevant Higgs-mass eigenstate has dominant components in the Higgs 
doublet $2 1^,1^ or, equivalently, when the matrix elements 02i and s^Oij related to the 



Higgs fields H(^4-i) are small. In the CP-conserving case, the Higgs boson that couples 



predominantly to the gauge bosons always fulfils these properties |^, i.e. O12 and O13 are 
zero for the CP-even Higgs states in this case. 

The general analytic expression for the effective Higgs-boson-bottom-quark coupling 



4T| may be obtained by considering the vertex graphs shown in Fig. ^. However, the mag- 
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Figure 20: Feynman graphs mediated by the exchange of (a) gluinos g and (b) Higgsinos 
hi2 that give rise to an effective one-loop $2^6 couphng. 



nitude of these corrections and the phases involved in these are strongly model-dependent. 
Furthermore, the gluino-exchange diagram in Fig. ^(a) usually represents the dominant 
contribution, which has a counterpart in the effective Higgs potential only at the two-loop 
level. To be precise, the effective one-loop Yukawa coupling of the h quark to the neutral 
field component of the $2 doublet, $2 = <^2 + ^^25 is given by 



C^oT.y. = Ahb^lbibR + H.c. 



-*066 



with 



Ah, 



2a. 



mgij,l{ml^,ml^,\mg\ 



+ 



h' 



jA/i/(m|,m|,|/in, 



(5.15) 



(5.16) 



where as = g'^/^in) is the SU(3)c fine structure constant, ht is the top Yukawa coupling, 
61 ("^ti) ^^^ ^2 ("^62) ^^^ ^^^ mass eigenvalues of the scalar-bottom (top) quarks, and 



mr 



I {a, b, c) is the one- loop function 
/(a, b, c) 



abln{a/b) + bcln{b/c) + acln(c/a) 



[a — b){b — c)[a — c) 

Note that the total bottom-mass corrections depend on the relative size of the gluino- 
mediated graphs to the Higgsino-mediated ones, shown in Figs. pO|(a) and ^(b), respec- 
tively. Therefore, as was mentioned above, the two different quantum corrections strongly 
depend on the model under study. In fact, it is straightforward to introduce these quan- 
tum effects in the analysis of the Higgs couplings. Taking both CP-violating vertex and 
Higgs-mixing effects into consideration, the effective Lagrangian for the Hibb couplings 
reads 



*-Hhb 



Y.H(i-i) b[[hb02i + Re(A/ib)03i - \ui{Ahb)cpOii 



i=l 



-i\^hi,Sf30ii - Re{Ahh)ci30ii - lm{Ahb)03i 75} &, 



(5.18) 
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where 






and Ahb/hb is given in Eq. (|5.16|) . 

As we have detailed above, the actual size of the CP- violating vertex effects depends 
both on 0CP and aig{mg^). The latter phase does not directly enter the calculation of 
the one-loop effective Higgs potential. To avoid excessive complication in the analysis, we 
have decided to present the results assuming that the vertex effects are very small, which 
is typically true for values of tan/3 < 4, like the ones considered here. For larger values of 
tan/5, instead, these effects can no longer be ignored. However, a specific model is then 
needed in order to be able to determine their significance. For example, for an appropriate 
choice of SUSY-mass parameters and phases, the vertex effects can even be tuned to zero 
independently of the CP-violating Higgs-mixing effects. 

6 Conclusions 

We have performed a systematic study of the Higgs sector of the MSSM with explicit 
CP violation, and analyzed the main phenomenological implications of such a theory for 
direct searches of Higgs bosons at LEP2 and the upgraded Tevatron collider. In this 
general theoretical framework, the tree-level CP invariance of the MSSM Higgs potential is 
considered to be sizeably broken by loop graphs involving trilinear CP-violating couplings 
of the scalar top and bottom quarks to Higgs bosons . These loop effects are taken into 
account by calculating the CP-violating RG-improved effective potential up to the next-to- 
leading order, in which two-loop leading logarithms of t-, 6- Yukawa and QCD corrections 
have been included 0. 

The analysis shows that the upper bound on the lightest Higgs-boson mass Mh^ 
obtained in the CP-violating MSSM is almost identical to that already derived in the CP- 
invariant theory, for both low and high tan/5 values. Nevertheless, the Higgs-sector CP 
non-conservation may drastically modify the couplings of the Hi scalar to the Z and W 
bosons. In fact, the production rate of the Hi boson is considerably affected at LEP2, 
for relatively light charged Higgs-boson masses, i.e. for 120 < Mh+ < 200 GeV, and for 
small and intermediate values of tan/3, i.e. 2 < tan/5 < 5. For larger values of tan/?, the 
soft-CP- violating couplings of the third generation are severely constrained by the two-loop 
SUSY Barr-Zee-type contribution to the electron EDM [^, leading to much weaker CP- 
odd effects. Because of the drastic modification of the HiZZ and H1H2Z couplings for 
low- and intermediate-tan /3 scenarios, we find that the current experimental lower bound 
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on the mass of the Hi particle may be dramatically relaxed up to the 60-GeV level in the 
presence of large CP violation in the Higgs sector of the MSSM (see Fig. ^(b)). Therefore, 
a combined experimental analysis is required by considering all possible reactions that 
contain the two lightest neutral Higgs bosons, Hi and H2, and the charged Higgs bosons 
H^ in the final state, which are produced either singly or in pairs. In this respect, the 
upgraded Tevatron collider and, to a greater extent, the LHC are the almost ideal places 
to explore more efficiently the parameter space of the CP-violating MSSM. 

Another consequence of such a minimal SUSY scenario of explicit CP violation is 
that the mass splitting between the two heaviest Higgs bosons H2 and H^ may be of order 
20% for heavy charged Higgs-boson masses, Mh± ^ 300 GeV (see also Fig. 0(b)). This 
is in agreement with earlier results reported in [^. Furthermore, we find that the strong 
mixing of the Higgs bosons may induce large CP violation in the vertices HiH^W^ , Hidd, 
H2dd, etc., that could even be of order unity. Most interestingly, CP violation can be 
resonantly enhanced in high-energy reactions mediated by scalar-pseudoscalar transitions 
that involve the nearly degenerate H2 and H^ states, especially when the mass difference 
Mff^ — Mh^ is comparable to the decay widths Th^ and Th^ |]T0[. Such resonant effects of 



CP violation may be tested at future high-energy colliders such as the LHC, NLC and the 
FMC 1 1^,11, 12, 13, 1^,15]. Finally, the analysis presented in this paper clearly demonstrates 



that the MSSM with explicit CP violation in the Higgs sector constitutes an interesting 
theoretical framework, which will have a significant impact on 5-meson decays, dark-matter 
searches and electroweak baryogenesis. 
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Note added 

While finalizing our paper, we became aware of two very recent works 0,5O| that treat 



some of the topics we have been studying here. In [^, the author mainly concentrated on 
the effect of CP-violating Higgs mixing on the Higgs-fermion-fermion coupling for relatively 
low values of /i, i.e. fi = 250 GeV, and Ma = 200 GeV. In this regime, CP violation in 
the sector of the lightest Higgs boson Hi is small. Instead, we effectively consider smaller 
values of Ma, i.e. 120 < Mjj+ < 180 GeV, and higher values of /i, which can lead to 
significantly larger CP-violating effects in the Hi sector. In this respect, our conclusions 
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differ from those presented in H3]; otlierwise, we find agreement with the results regarding 



the -^2// couphng. The authors of Ref. [Q discuss Higgs-boson production cross sections 
for a future e^e~ NLC in a general CP-violating two-Higgs-doublet model. Here, instead, 
we are mainly interested in possible effects at LEP2 and the upgraded Tevatron collider, 
within the MSSM with radiatively induced CP violation in the Higgs sector. In addition 
to [^9| , p0| , we present analytic expressions for the Higgs-boson masses and mixing angles, 
and pay particular attention to the EDM constraints. 
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A Analytic expressions of quartic couplings 

The dominant contribution of radiative interactions to quartic couplings comes from en- 



hanced Yukawa couphngs of the third generation. The relevant Lagrangians |^ , including 
CP-violating sources, are given by 



-C 
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(A.2) 

(A.3) 
(A.4) 



with Q^ = (tL,bL), U* = iji, D* = h^. In addition, we also consider next-to-leading order 



QCD quantum corrections. These corrections have been computed in full detail in |17,18 



and may easily be implemented in the analysis, as they form a CP-invariant subset of 
graphs by themselves. 

Employing the interaction Lagrangians ( |A.l|) -( |Al4|) , it is straightforward to extend 
the two-loop analytic results in [^ to the case of CP violation. In this way, we find 
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In Eq. ([A.12| ), rrit is the top-quark pole mass, which is related to the on-shell running mass 
rrit through 

mt{mt) = , \_ . (A.14) 

It is important to remark that the above expressions are not equivalent to the ones that 
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would be obtained by taking the expressions given in Ref. [Q, and considering all mixing 
parameters to be complex. If this were done, incorrect results would be obtained. 

The RG analysis under consideration assumes a single step decoupling of the scalar- 
quark fields. This assumption is only valid if the mass splitting among the scalar-quark mass 
eigenstates is relatively small 0. Specifically, the expansion becomes more trustworthy if 

2 2 

mf — mf 

-^ ^ < 0.5, (A.15) 

m- + mj 

where m? and m? are the squared mass eigenvalues of the scalar-top quarks. This last 
restriction also applies to other Higgs-mass analyses that have been performed at the next- 
to- leading order @J^,^ . If one assumes that the approximate inequality (|A.15|) holds true. 



the scale Mg^gy "^^J then be safely defined as the arithmetic average of the scalar-top mass 
eigenvalues squared. 



-'^■'SUSY 2 



Uml+ml). (A.16) 



One should bear in mind that our RG analysis also relies on an expansion of the effective 
Higgs potential up to operators of dimension 4. The contribution of higher- dimensional 
operators may only be neglected if 2\mtAt\ ^ Ml^gy ^^^ 2|mt/i| < tan/3M|ugY- Moreover, 
in the high tan/3 regime, in which the bottom- Yukawa interactions become more relevant, 
we have assumed that the scalar-bottom masses are of the order of the scalar-top ones and 
that similar bounds on their respective mixing mass parameters are fulfilled. Observe that, 
in order to evaluate Higgs-boson pole masses, Higgs-vacuum-polarization contributions 
should be included in the calculation. In general, for our choice of the renormalization 
scale, the vacuum-polarization contributions are small. These contributions would only 
be necessary in a calculation of Higgs-boson masses that goes beyond the approximation 
presented here. Therefore, the present approach is very analogous to that studied in Ref. 
|5^ in the MSSM framework with CP-invariant Higgs potential. 

The computation of the Higgs-boson mass-matrix elements considered here and in 
0] is still affected by theoretical uncertainties, most noticeably, those associated with two- 
loop, finite threshold corrections to the effective quartic couplings of the Higgs potential. 
Recently, a partial, diagrammatic, two-loop calculation of the mass of the lightest CP-even 



Higgs boson has been carried out |53|. In the limit of large Mh+, the additional two- 
loop threshold corrections lead to a slight modification of the dependence of the lightest 
CP-even Higgs-boson mass on the scalar-top mixing parameters. For instance, although 
the upper bound on the lightest CP-even Higgs-boson mass for TeV scalar-quark masses 
is approximately equal to that obtained through next-to-leading order, which is also the 
approach we followed in the present work, the upper bound on the Higgs-boson mass is 
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reached for slightly different values of \At\, i.e. \At\ ~ 2Msusy instead of \At\ = \/6Msusy, 
with a weak dependence on the sign of Af. Similar results were obtained by means of a 
two- loop calculation of the effective potential |S^ . Nevertheless, a complete diagrammatic 
analysis of the two-loop corrections induced by the t- Yukawa coupling, which are included 
at the leading-logarithmic level in our computation, is still lacking. 
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